Numerical studies of confinement in the lattice Landau gauge 
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Critical conjectures on confinement in the Landau gauge is numerically tested in focus to Gribov copy effects. 
One of the subjects is of the Kugo-Ojima confinement criterion and the other is of various viewpoints in the 
Gribov-Zwanziger theory. We use the smearing gauge as a reference gauge free of Gribov copy, and performed 
three types of simulations, log U, (7-linear and log U in the smearing gauge. It is found that Gribov copy effect 
on the Kugo-Ojima parameter is small, log U and [/-linear simulations yield only global scale factor difference 
in gluon propagator and in ghost propagator, and about 10% difference in Kugo-Ojima parameter. The horizon 
function defined by Zwanziger is evaluated in three types of gauge field and compared. All data show the negative 
horizon function as expected. 



1. INTRODUCTION 

Gribov pointed out that there exists a funda- 
mental problem with respect to the field theo- 
retic description of QCD in its quantization, due 
to gauge fixing degeneracy H. Gribov conjec- 
tured that the linear potential of quarks may be 
caused by an enhancement of the singularity of 
the ghost propagator due to the restiriction of the 
gauge field on the transverse plane. Zwanziger 
developed further the lattice gauge theory along 
the line of Gribov, and discussed its continuum 
limit g. 

Kugo ond Ojima developed the operator for- 
malism of QCD starting from the standard 
BRST (Bccchi-Rouet-Stora-Tyutin) symmetric 
Lagrangian, and derived Kugo-Ojima criterion, 
a sufficient condition for the colour confinement, 
i.e., the absence of free single coloured states in 
the asymptotic Hilbert space J^]. Even though 
there is no question in validity of the starting 
Lagrangian as the perturbation theory, the va- 
lidity is unclear as the full theory, if one thinks of 
the Gribov problem which emerges in the 'large' 



amplitude region of the gauge field. There is a 
discussion of Fujikawa || and Hirschfeld|| trying 
to justify the standard Lagrangian. This point 
is intimately related to a problem in numerical 
measurement [Q of gauge non-invariant quantity 
in presence of Gribov copyg. 

The Kugo-Ojima theory and Gribov-Zwanziger 
theory appear to stand on different footings, but 
each presents its own critical conjectures on con- 
finement problems. It is meaningful to explore 
their validity numerically by lattice QCD dynam- 
ics. We test how numerically satisfied is the 
Kugo-Ojima criterion that a coefficient u£(0) in 
the two-point function produced by the ghost, the 
antighost and the gauge field becomes —5%, and 
investigate how gluon propagator and ghost prop- 
agator behave in the infrared region. We also look 
at value of Zwanziger's horizon function to test if 
his horizon condition is satisfied ||. The theoreti- 
cal problem of quantization caused by the Gribov 
copy reflects as the Gribov noise problem in mea- 
surement of gauge non-invariant quantity, and it 
is favoured to choose a unique copy among oth- 
ers, i.e., a new gauge without Gribov copy 
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(see [|| in more detail). As finding algorithms of 
unique gauge fixing is difficult in various gauges in 
general, random statistical method is commonly 
used, but it is cautioned that its performance 
is sensitive to incorporated algorithms in some 
cases Laplacian gauge fixing is well known as 
a unique gauge fixing, but only known method 
in Landau gauge is smearing gauge fixing We 
call this gauge as smearing gauge and adopt it as 
a reference gauge to investigate the Gribov copy 
effect. 

There are some options of gauge field defini- 
tion A^U) in lattice gauge theory. We compare 
numerical data with respect to two types of defi- 
nition, log U and [/-linear. 

2. Gauge fields and lattice Landau gauge 

Options of A^(U) for SU(3) lattice QCD are 

1. logZJ; U x ^ = e A *-"-, = -A x ^, where 
eigenvalue of A x ^\ < 47r/3 [px| . 

2. [/-linear; A Xjfl = ±(U Xtli - C^, M )| tP i.@. 

In both cases the Landau gauge, dA 9 = 0, can 
be characterized in use of optimizing functions 
Fu{g) of g, such that SFjj{g) = for any Sg, 
respectively; (trl means the traceless part) 



1. F u (g) = \\A*\\ 2 = j: x , fl tr(As x J 



2- ^(9)=Ex,m( 1 



|Retr[/^) 



In both options, the variation of the optimizing 
function, Fu(g), under infinitesimal gauge trans- 
formation g _1 Sg = e, reads as 

AFu(g) = -2(dA9\e) + (e\ - dD{A*)\e) + • ■ ■ , (1) 

where D^A) denotes the covariant derivative in 
each definition respectively; in short notations, 
U = U x>fl , A = A x ^, and d^tj) = (f>(x + fi) - (f>{x), 



1. D ll {A)cl> = S{A)d IJt (j>+[A,cj>}, 
where A = adjA — [A, •], S(x) = 



2. D^U)ct> = - 



1 (U + W 



,9a<t> 



(2) 

x/2 
thO/2)' 

+[A0]-(3) 



3. The Kugo-Ojima confinement criterion 
and the Gribov-Zwanziger's theory 

3.1. The Kugo-Ojima criterion||2[ 

The Kugo-Ojima criterion is stated as follows; 
u ab (0) = -6 ab , where u ab is defined by the fol- 
lowing two-point function as, 



KV) 



V ^ 

x,y 



)(4) 



xy 



trl. 



where angle brackets denote sample average, and 
A a is a normalized (trA a ^A b = d ab ) antihermitian 
basis of Lie algebra, V a lattice volume. 

3.2. Zwanziger's theory^ 

Zwanziger classified transverse space of gauge 
field into some sets. Fjj(g) takes the local mini- 
mum in the Gribov region f2, and the global min- 
imum in the fundamental modular region A, 
A C fi. Within A, there exists no gauge fix- 
ing degeneracy up to global gauge transforma- 
tion. We call this gauge A gauge. For purpose of 
gauging away Bloch wave states, Zwanziger nar- 
rowed down A to core region 5 C A : Let the 
space of configuration of period L be 11^, and 
the fundamental modular region in the Hl he 
Al, and let A^ = Anl n Hl, then the core re- 
gion is given by = A£\ Note, however, that 
by this definition, considerably large amount of 
points in Al fail to come into A^ not because 
they can be gauge transformed to the point in 
A^ , but because they fail to satisfy L-periodicity 
condition when gauge transformed to Ajvl- From 
this view point, narrowing down to the core re- 
gion looks like a kind of artificial selection of 
configurations. Thus Zwanziger's basic hypoth- 
esis is that partition functions defined by path 
integrals over and A/,, give the same limit, 
liniL^oo Ze l = limL_,oo Z\ L . Since this hypothe- 
sis is highly dynamical, it is to be checked some- 
how, e.g., by simulation. One can show that the 
following horizon function H(U) is negative in 
the above defined core region S by estimating in- 
crease of Fjj(g) by Bloch wave excitation. H(U) 
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is given as follows; 

H(u) = J2 G ^y aa - * E W> 



x,y,a 



Gfj,uxy = tr I A" T 25„— (-AOA 1 
with 

1. E(U) = g£> (\ at S(Ai)\ a 

l,a 

2. £(£/) = J2 ^ trC/i, 



(5) 
(6) 

(7) 
(8) 



Statistical average defines a tensor G^ vxy as 
(Gfflxy) — G^ vxy 8 ab , provided colour symmetry 
is not broken. The Fourier transform of the ten- 
sor G^xy takes a form 



G 



where e = (E(U))/V. It is related with the hori- 
zon function as 



(H(U)) 
V 



lim G MM (p) 



(10) 



Note that e = 4 if all link Ui = 1. For pur- 
pose of path integral estimation of the parti- 
tion function, a set augmented core region 
* = {[/: H(U) < 0} n fi ( 3 C * C n ) is de- 
fined, and from estimation of Z*, one derives the 



infinite volume limit Q that limy^oo — — — = 

which is called horizon condition. Putting Kugo- 
Ojima parameter as u ab (0) = —6 c, one finds 
from (||), (|ic|), that the horizon condition is writ- 
ten as (j^J + 3c — e = 3 — — ^ =0. 



4. Method of Landau gauge fixing, numer- 
ical results and discussion 

Our simulation was done for (3 = 6.0 in lat- 
tice size 16 4 . Simulation of gauge fixed theory 
can be simply executed by gauge transforming a 
Boltzmann sample to the gauge, and measuring 
a quantity in question 0. This algorithm is un- 
ambiguous and correct only if the gauge is free 



of Gribov copy, or the quantity is gauge invari- 
ant. It is easily understood that if the quantity 
is gauge non-invariant and if there exists Gribov 
copy, numerical results depend not only on the 
gauge itself, but also on the algorithm employed 
which brings a sample to a gauge copy||. We 
performed three types of Landau gauge calcula- 
tions, combinations of gauge field definitions and 
gauge fixing algorithms. 

1. Gauge fixing with \ogU option using a 
direct Newton method, i.e., gauge trans- 
formation phase 6 being given by <fi = 
(-dD(A))~ 1 dA^^. Gribov copy exists. 

2. Gauge fixing with logt/ option using the 
smearing gauge fixing jj]. This smearing 
gauge is a unique Landau gauge. But our 
execution shows that the smearing gauge at 
j3 = 6, 16 4 does not coincide completely 
with A gauge. 

3. Gauge fixing with [/-linear option using 
the non-stochastic-over-relaxation method 
after the process 1. Gribov copy exists. 

In our smearing gauge fixing, accuracy of pla- 
quette value E p of smeared vacuum configuration 
is set as 1 — E p < 10~ 7 , since there appear a long 
plateau around 1—E p ~ 10~ 4 in smearing process 
for many samples, and the finally obtained gauge 
fixed configuration depends critically on the sit- 
uation if the smearing is finished before or after 
the long plateau. 

Kugo-Ojima parameter is found to be uf(0) ~ 
—0.65 for cases 1,2. and there is no remarkable 
difference. The magnitude of the parameter is 
about 10% less in case 3. u% does not look like 
tending to -1 up to the lattice size 16 4 . The reason 
is unclear so far. And Table 1 shows the Kugo- 
Ojima parameters and values related to the hori- 
zon function. All values indicate negative horizon 
function|]l2||. 

Gluon propagator is infrared finite, and 
log — log plot of unrenormalized gluon propaga- 
tor is given in Fig 1, and about 20% difference 
of case 3 appear as a global scale factor, i.e., 
wave function renormalization, and this situation 
agrees with the discussion in the literature flli"||. 
Case 2 is almost same as case 1. 
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Table 1 

Results of the Kugo-Ojima parameter c, 'trace' e divided by the dimension 4, and h = c — e/4. 



gauge fixing 


c 


e/4 


h 


log U 


0.628(94) 


0.943(1) 


-0.32 


smeared log U 


0.647(101) 


0.943(1) 


-0.30 


[/-linear 


0.576(79) 


0.860(1) 


-0.28 
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Figure 1. The log of the unrenormalized gluon 
propagator multiplied by (qa) 2 as a function of 
the log qa where q is the lattice momentum. Tri- 
angles and the fitted line —0.743 log qa + 0.576 
are in the log U , and diamonds and the fitted line 
—0.747 log qa + 0.382 are in the [/-linear. 



Figure 2. The log of the unrenormalized ghost 
propagator multiplied by (qa) 2 as a function of 
the log qa where q is the lattice momentum. Tri- 
angles and diamonds are the same as Fig[l]. The 
fitted line is -0.390 log ga + 0.278 in the log U and 
-0.373 log qa + 0.380 in the [/-linear. 



log — log plot of ghost propagator is given in Fig 
2. The difference of case 1 and case 3 is the global 
scale factor as in the case of gluon propagator. 
Case 2 is almost same as case 1. 
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